Scaling of electrical and thermal conductivities in an almost integrable chain 
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Many low-dimensional materials are well described by integrable one-dimensional models such as 
the Hubbard model of electrons or the Heisenberg model of spins. However, the small perturbations 
to these models required to describe real materials are expected to have singular effects on transport 
quantities: integrable models often support dissipationless transport, while weak non-integrable 
terms lead to finite conductivities. We use matrix-product-state methods to obtain quantitative 
values of spin/electrical and thermal conductivities in an almost integrable gap less chain (an XXZ 
spin chain with staggered fields, or equivalently a spinless fermion chain with staggered on-site 
potentials). The results at low temperatures validate a scaling theory based on bosonization. 



The physics of many one-dimensional systems with ide- 
alized interactions is rather special: the quantum Hamil- 
tonian has infinitely many independent conserved quanti- 
ties that are sums of local operators. Such Hamiltonians 
are called "integrable" in analogy with classical Hamil- 
tonian systems that decompose into independent action 
and angle variables. Examples relevant to experiments on 
crystalline materials include the Hubbard model of elec- 
trons and the XXZ model of spins; ultracold atomic sys- 
tems can realize integrable continuum models of bosons. 
However, in all these cases it is expected that integrabil- 
ity is only an approximation to reality and that experi- 
mental systems have integrability-breaking perturbations 
which, while small, drastically change some of the phys- 
ical properties. 

Transport of a conserved quantity such as charge, spin, 
or energy provides an experimentally important example 
of the effects of non-integrable perturbations. Integrable 
systems often have dissipationless transport or a "Drude 
weight" D in their frequency-dependent conductivity: [l|- 

H 



a(oj) = 2itD6(uj) 



(1) 



where the ellipsis represents the nonsingular part of the 
conductivity. An experimental example of similar physics 
is the large, anisotropic thermal transport observed in 
Sri4Cu2404i and attributed to a long mean free path 
of quasi- ID magnons [16 18]. We are interested here in 
the question of how the presence of generic integrability- 
breaking terms regularizes the zero-frequency conductiv- 
ity 0-0, 12, 14, 1§-23|, as we expect many quasi-one- 
dimensional systems to be well described by integrable 
Hamiltonians plus weak non-integrable perturbations. 

We apply a combination of numerical (matrix prod- 
uct states) methods to a suitable gapless non-integrable 
Hamiltonian to observe quantitatively the destruction of 
the thermal and electrical Drude weights. We compute 
the corresponding conductivities and analyze how they 
depend on temperature and the strength of the non- 
integrable perturbation. The results are compared to a 
scaling theory based on bosonization. One might expect 



that on adding non-integrable terms, which generically 
destroy conserved quantities and hence the Drude weight, 
bosonization allows to obtain a parameter-free prediction 
for the conductivity at low temperature. By comparing 
with our numerics we obtain an understanding of the 
accuracy and range of validity of the bosonization ap- 
proach. 

We employ the XXZ chain in presence of a staggered 
magnetic field h which breaks integrability and focus on 
the parameter regime when this perturbation is irrele- 
vant in the scaling sense and does not open a gap. We 
expect our key results (such as the scaling behavior of the 
conductivity at low temperatures) to hold for any gap- 
less non-integrable model in which no conserved quantity 
protects the current. In the regime where h is relevant, 
its effects have been studied perturbatively [13, Hit an< l 
compared to experiments on spin diffusion in copper ben- 
zoate 0, H(| . A classic probe of quantum integrability, 
the absence of level repulsion in the energy spectrum, is 
used to show that this model is indeed non-integrable. 
We obtain numerical results for the electrical/spin and 
thermal conductivities using density matrix renormal- 
ization group (DMRG) methods which were developed 
in the past few years to access finite-temperature dy- 
namics of correlated systems. Our model involves only 
nearest-neighbor interactions and can thus be studied ef- 
ficiently via DMRG. By adapting previous bosonization 
techniques [H, [l| to the non-integrable perturbation in 



0.5 




/, (0.99,0.4^5) 


0.4 


AF 


(0.95,0.4^/ 


0.3 




/ ■ (0.9,0.31) / 


0.2 




/ / / 
' -(0.85,0.2) / 


0.1 




r LL / FE 









0.9 

-A 



FIG. 1. The phase diagram [H, El of © . The points share 
the same Luttinger liquid parameter K « 2.4 computed by 
DMRG [H]. 



2 



500 
400 
300 
200 
100 


i 

500 
400 
300 
200 
100 



Lh=0.0 




FIG. 2. The distributions of {E l+1 - Ei}, where E x < E 2 < 
. . . < -B4862 are the eigenenergies of (J2j for L = 18, A = —0.8 
in the S z = 1, k = 2n/9 sector (the total magnetization S" = 
XL = i Sf an< i the momentum k are conserved). The green 
curve and red curve are best-fit exponential and Wigner- 
Dyson (orthogonal ensemble) distributions respectively. Nei- 
ther the Poisson nor the Wigner-Dyson distribution appears 
clearly if we do not restrict to a symmetry sector of the model. 
The crossover to Wigner-Dyson with increasing h is observed 
independent of the choice of A and symmetry sector. 



our model, we find testable predictions for the electrical 
conductivity at low temperature and compare them to 
the numerics. 

The model. The XXZ model in a staggered magnetic 
field is given by H = hi with 

hi = S*S* +1 + SfSf +1 + ASiSi +1 + (-If hSl. (2) 

Without the staggered field, ([2]) is gapless for | A| < 1 and 
gapped for | A| > 1. Bosonization leads to the low-energy 
effective Hamiltonian for |A| < 1 and infinitesimal h: 

H = i / dx ( n2 + ^ dx< ^ 2 ) + ch J dxcos (^^k^ 

H~ ^umklapp -*?band curv. -"higher terms in h ) 

(3) 

where II is the conjugate momentum of the bosonic field 
(f> with the canonical commutation relation [<p(x),H(y)] = 
iS(x — y). The first term in © describes a Luttinger liq- 
uid. The Luttinger liquid parameter K is given through 
Bethe ansatz: 2K arccos(— A) = ir, and the coefficients 



v and c are also known exactly [31], l32j . As the scaling 
dimension of h is 2 — K, the second term in ([3]) is relevant 
and opens a gap for K < 2 or — V2/2 < A < 1; the term 
is irrelevant and ([2]) remains in the gapless Luttinger liq- 
uid phase for K > 2 or -1< A < -V2/2 (Fig. [TJ. 

Integrability is well-defined in classical mechanics, but 
the definition of its quantum counterpart remains a sub- 
ject of debate [33| . It is generally believed that the level 
spacing distribution (the distribution of the differences 
of the adjacent eigenenergies) is the exponential distribu- 
tion for an integrable model, as levels appear as a Poisson 



process, and the Wigner-Dyson distribution for a nonin- 
tcgrable model. Intuitively, two nearby levels in an inte- 
grable model likely have different values of at least one 
integrable quantity, and thus live in different sectors of 
Hilbert space that are independent of each other; hence 
their energies are uncorrelated. Non-integrable models 
do not have an extensive number of such sectors and 
show energy level repulsion. The belief has been verified 
numerically on a variety of models [13, EH . We perform 
an exact diagonalization of (f2|) with periodic boundary 
conditions. Fig. [5] shows the level spacing distributions, 
and the crossover from Poissonian behavior at h = 
to the Wigner-Dyson distribution at nonzero h is clear. 
Hence ^ is nonintegrable for nonzero h. 

Numerical approach. The DC charge (c) and heat (h) 
conductivities can be computed via the Kubo formula 

a = lim lim -!- Re / (J(t)J(Q)) dt , (4) 

tjkf— too L— >oo LT Jq 

where the corresponding current operators J = Yl%=iJi 
are defined through a continuity equation 



i=2 
L 



(5) 



dtS- = j c .i - ici+i => J c = i^2[hi-i,S*] . 



i=2 



The current correlation functions (J(t)J(O)) can be com- 
puted efficiently u sing the real-time finite-temperature 
DMRG algorithm [3^ introduced in [ll|. DMRG is es- 
sentially controlled by the so-called discarded weight e 
(we ensure that e is chosen small enough and that L 
is chosen large enough to obtain numerically-exact re- 
sults in the thermodynamic limit). The simulation is 
stopped when the DMRG block Hilbert space dimension 
X reaches x ~ 1000 — 1500. This allows to access time 
scales t ~ 10—20 which are much larger than the inherent 
microscopic scale t = 1. 

Results for (J(t)J(0)) are shown in Fig. [3] In the in- 
tegrable case h = 0, the heat and charge Drude weights 



D 



lim lim 

t— >co L— >oc 



Re (J(t)J(O)) 
2LT 



(0) 



are finite. A nonzero h > renders the model nonin- 
tegrable; the current correlation functions decay to zero 
at large times, and the conductivities become finite. In 
order to compute the integral in Eq. (j4|) quantitatively, 
(J(t)J(O)) needs to be extrapolated. The heat current 
correlation function at intermediate to large tempera- 
tures 0.5 < T < oo (see Fig. [3jb)) can be fitted by a 
single exponential function exp(— At) j3?| • Oscillations 
develop at small T, but it is reasonable to assume that 
(J(t)J(O)) can be described by sums of exponentially de- 
caying terms exp(— \ n t + iu) n t) (the same holds for the 
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FIG. 3. Heat and charge current correlation functions 
( J(t) J(0)} at A = —0.85. Their Fourier transform determines 
the corresponding conductivities through Eq. ((4]). The Drude 
weight (Eq. JSJ) is nonzero only in the integrable case h = 0. 
A nonintegrable perturbation h > renders the conductivity 
finite. DMRG data (solid lines) is extrapolated using linear 
prediction (dashed lines). 

charge current correlation function). This motivates us 
to use so-called linear prediction |36| as an extrapolation 
procedure. Its stability can be tested by varying fit pa- 
rameters (e.g. the number of terms or the fit interval) 
and by checking sum rules (see below); we can obtain 
accurate results for the heat conductivity at any h and 
temperatures 0.2 < T < oo as well as for the charge con- 
ductivity at intermediate to large h and small T < 0.3. 

Scaling of the conductivities: numerical results. 
DMRG data for the heat and charge conductivities is 
shown in Fig. ffl For fixed T and small h, one expects a 
to diverge as [23| 

a ~ h,- 2 for h -> . (7) 

This is confirmed by our results (for the thermal case see 
the inset to Fig. [Ha)). At small T, a c features power 
laws with nontrivial exponents: 

cr c - T a " for T -> . (8) 

Our model is a Luttinger liquid at low energies and one 
thus expects a c to be a universal function of the Luttinger 
liquid parameter K only. This is verified in Fig. E^b) 
which shows cr c (T) for different parameter sets (A, h) 
which share the same K (K can be determined from an 
independent ground state DMRG calculation [111). The 
exponent a c (K) varies strongly with K; it is consistent 
with the analytic predition a c = 3 — 2K established be- 
low (see the insets to Fig. BJb)). The heat conductivity 



is only accessible for intermediate temperatures T > 0.2; 
our data in this regime is almost independent of K (see 
Fig.dja)), suggesting that we have not reached the limit 
of low T. 

Since for h = the heat current is conserved by 
the Hamiltonian, the AC conductivity o~h(tjJ,h = 0) = 
2irD5(ui) features a Drude peak only. By generalizing 
Eq. dl]) to finite frequencies d, 0, we can straightfor- 
wardly compute 07, (w, h) and demonstrate that it indeed 
becomes a (^-function series for h — > 0; the frequency- 
integrated heat conductivity just yields the Drude weight 
('sum rule'). This is illustrated in Fig. [5]and provides an 
independent test for the reliability of our extrapolation 
procedure. 

Scaling of the charge conductivity: theoretical analy- 
sis. We now present an analytic calculation of the low- 
temperature behavior of the charge conductivity a c as- 
suming that the terms in Eq. ([3]) generated by h are al- 
ready sufficient to destroy integrability. While this de- 
scribes our data well, another possibility [2(| is that in- 
tegrability breaking results from the combination of h 
and umklapp, which leads to non-commuting h — > and 
T — > limits; this would also explain our data if we are 
currently in the small- h limit, but would predict a differ- 
ent behavior to set in at even lower temperatures. We 
start with a simple scaling analysis. Combining several 
simple assumptions, which are likely to hold for other 
ID models, we establish a scaling form for the conduc- 
tivity in which all the exponents are determined up to 
one number, the scaling dimension of the integrability- 
brcaking perturbation. For our model, bosonization will 
verify this scaling ansatz and yield the scaling dimension. 

Scaling analysis. It is reasonable to assume that 
Re(J(t)J(0))/LT w A(A,h,T)exp(-j(A,h,T)t) at 
long time as correlations typically decay exponentially at 
finite T. The oscillation of this correlation function (see 
Fig. [3]) of Re( J(t) J(0)) /LT is not taken into account, as 
it cancels out in computing the integral We also as- 
sume the amplitude A(A, h, T) does not vanish as T — > 
(note that the Drude weight D(A,h = 0, T) is nonzero 
and continuous as T — > 0). Then, (H} implies er c ~ 7" 1 , 
and o~ c takes the scaling form 

a c {A, h, T) = f{A/T^\h/T^)/T, (9) 

where [A] and [h] are the scaling dimensions of A and h, 
respectively. Note that er c ~ T _1 for [A] = [h] = or at 
the phase transition A = — v2/2, h = o(l). 

In the perturbativc regime (i.e., infinitesimal h), 
[A] = as there is no renormalization of A (it is ex- 
actly marginal). Then ([9]) simplifies to a c (A,h,T) = 
f(A,h/T [h] )/T. As one expects a to diverge as h~ 2 by 
a golden-rule argument j23J unless this perturbation is 
inefficient in inducing scattering, we take f(x) ~ x~ 2 . 
Then 

g c ~ h~ 2 T 2 W- 1 = h- 2 T 3 ~ 2K (10) 
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FIG. 4. Scaling of the heat and charge conductivities of a nonintegrable gapless spin chain. At low temperatures, a c (T) ~ j^ei-Kj 
is governed by power laws whose exponents are a universal function of the Luttinger liquid parameter K only. At fixed T, a 
diverges as hT 2 for h — > with h being the strength of the nonintegrable perturbation j23|]. Note that our definition of an via 
Eq. Q differs from the most common one [TJ by a factor T -1 . 



where in the second equality we have substituted the 
bosonization result [h] = 2 — K. Note that a c ~ T^ 1 for 
K = 2 or at the phase transition A = —y/2/2, h = o(l), 
consistent with the result at the end of the last para- 
graph. 

The results of scaling for infinitesimal h will be repro- 
duced by a perturbative field theory calculation using 
bosonization, which more convincingly justifies the as- 
sumptions in this section. However, it is worth pointing 
out that from the scaling analysis we still expect scaling 
of conductivity at low temperature in a gapless ID sys- 
tem even when bosonization is inapplicable. In general 
a gapless ID system with a single velocity of low-energy 
excitations will be described by a conformal field theory 
(CFT) at long distances, and such theories are effectively 
ballistic as left-moving and right-moving excitations de- 
couple. We expect that the basic picture that conductiv- 
ity is controlled by the leading irrelevant operator that 
induces scattering will still apply even when the CFT is 
more complicated than a free boson. 

Bosonization. We now sketch a bosonization approach 
closely following Refs. H and H which studied transport 
in the integrable XXZ model. Note that in our case 
no conserved quantity protects the Drude weight. The 
current operator (O reads J c = —vy/K/n J dxll, and 
the Kubo formula for the AC conductivity reduces to 
a c (u>) = iK w ((fx/)) r (u>)/n. We perform a perturbative 
field theory calculation of the retarded correlation func- 
tion (4><j)) r to leading order in h, in -ff U mkiapp, and in 
-ffband curv. • The leading term governing the DC conduc- 
tivity reads a c = h~ 2 T 3 ~ 2K /C(K), which is consistent 
with Eq. (JTUJ) j38|. This term can be attributed to the 
staggered field, which is reasonable since the umklapp 
term alone does not destroy integrability and is thus less 
effective in scattering currents. 

Outlook. Our work demonstrates an approach valid 



for many actual ID materials, in which integrability- 
breaking terms are likely to be present but small. We 
studied one specific model but expect that our key result 
- a power-law scaling of the conductivity a ~ T a with 
a universal exponent determined by the Luttinger liquid 
parameter - should be a general qualitative feature of 
any gapless nonintegrable model in which no conserved 
operator protects the current. Quantitative results for 
other nonintegrable perturbations can be obtained by 
the numerical framework used in this paper. It should 
be possible to compute optical charge conductivities for 
comparison to experiments on conducting polymers and 
other systems. 

On a more basic level, quantum critical transport in 
one dimension is controlled by the leading irrelevant op- 
erators if and only if those destroy integrability. In higher 
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FIG. 5. AC heat conductivity. A Drude peak emerges as 
h —¥ 0: The integrated conductivity is independent of h (right 
inset) and equal to the Drude weight (left inset). 
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dimensions, quantum critical transport is different be- 
cause the critical theory is believed to be non-intcgrablc, 
and transport properties are actively being studied by 
methods from high-energy physics. Our results provide 
a constraint on these methods in a case where direct com- 
putation of transport coefficients is possible. 
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